Nuclear-spin relaxation of ^°^Pb in ferroelectric powders 
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Motivated by a recent proposal by O. P. Sushkov and co-workers to search for a P,T-violating 
SchifT moment of the ^°^Pb nucleus in a ferroelectric solid, we have carried out a high-field nuclear 
magnetic resonance study of the longitudinal and transverse spin relaxation of the lead nuclei from 
room temperature down to 10 K for powder samples of lead titanate (PT), lead zirconium titanate 
(PZT), and a PT monocrystal. For all powder samples and independently of temperature, transverse 
relaxation times were found to be T2 « 1.5 ms, while the longitudinal relaxation times exhibited a 
temperature dependence, with Ti of over an hour at the lowest temperatures, decreasing to Ti ~ 7 s 
at room temperature. At high temperatures, the observed behavior is consistent with a two-phonon 
Raman process, while in the low temperature limit, the relaxation appears to be dominated by a 
single-phonon (direct) process involving magnetic impurities. This is the first study of temperature- 
dependent nuclear-spin relaxation in PT and PZT ferroelectrics at such low temperatures. We 
discuss the implications of the results for the Schiff-moment search. 



PACS numbers: ll.30.Er, 67.80. Jd 

I. INTRODUCTION 

The ^°^Pb nuclear system (nuclear spin I — 1/2; mag- 
netic moment /i « 0.58 hm] isotopic abundance ~ 22%) 
in ferroelectric solids has been proposed for a search for a 
SchifF moment associated with simultaneous violation of 
parity (P) and time-reversal invariance (T) in fundamen- 
tal interactions [T] (see also a discussion of the sensitivity 
of such search in Ref. ^). The idea is that, due to the 
SchifF moment, a ferroelectric sample would acquire a 
P,T-odd magnetic polarization along the direction of its 
electric polarization. 

In a compound free of unpaired electrons, at low 
enough temperatures when lattice vibrations freeze out, 
we have an isolated system of nuclear spins (in the sense 
that the spins can freely rotate, but the nuclei cannot 
move with respect to the lattice) that interact among 
themselves solely via magnetic dipole-dipole interactions. 
It was shown in Ref. [T] that lead nuclei in perovskite 
ferroelectric compounds such as lead titanate (PT) ex- 
perience large effective electric field, producing P,T-odd 
energy shift. The effect is enhanced in ferroelectric lead 
compounds due, in part, to the rapid (faster than Z"^) 
scaling of the effective electric field with the charge of 
the nucleus, Z, which is also combined with the ~ Z^^^ 
scaling of the Schiff moment, see, e.g., Ref. [3, Sect. 
8.5.1]. Thus, the ferroelectric compounds proposed for 
the Schiff-moment search combine the advantage of a 
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large number of particles accessible to a solid-state ex- 
periment with an enhancement of the effective electric 
field similar to that in diatomic molecules (see, for exam- 
ple, the review [31). 

A crucial requirement for the success of the P,T- 
violation experiment is a detailed understanding of the 
behavior of the isolated nuclear-spin system, and in par- 
ticular, its relaxation properties. The connection of 
the relaxation and fundamental limits of sensitivity of 
condensed-matter P- and T-violation experiments was 
discussed in Ref. [2] , and a related question of magnetic 
susceptibility was discussed theoretically and modeled 
numerically in Ref. [5j. 

In this work, we begin experimental investigation of 
the lead nuclear-spin system in the context of this ap- 
plication by studying, using nuclear magnetic resonance 
(NMR) techniques, longitudinal (Ti) and transverse (T2) 
relaxation in ferroelectric powders and a PT monocrystal 
as a function of temperature in the range 10-293 K. We 
find that, at high temperatures, these ferroelectrics un- 
dergo magnetic relaxation due to a two-phonon Raman 
process. At temperatures much lower than the Debye 
temperature, relaxation is dominated by impurities and a 
single-phonon (direct) process. The obtained results en- 
courage us to further pursue experimental work towards 
a Schiff-moment search in lead-based solid ferroelectrics. 



II. PREVIOUS WORK 
A. Studies of lead NMR 

The body of work on solid-state NMR of lead prior to 
2002 was reviewed in Ref. [6] . Because of the large atomic 
number of lead, its compounds often display large chemi- 
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cal shifts, up to several thousand ppm. The isotropic and 
anisotropic parts of the chemical shift tensor, as well as 
their temperature dependences have been shown in Ref. 
[7J to be sensitive probes of the crystalline structure. The 
207pb sohd-state NMR study of PT [7^ was conducted at 
sample temperatures between -150 and 60°C. In Ref. [8], 
Pb NMR of a single-crystal ferroelectric lead magnesium 
niobate (PMN) was studied at temperatures down to 15 
K. The authors investigated the effect of electric poling 
on the NMR spectrum. They observed longitudinal re- 
laxation times for lead in the range 1-10 s. 

In some cases, NMR spectroscopy of ^"^Pb can be 
done in conjunction with the spectroscopy of other 
nonzero-spin nuclei in the same compound. For exam- 
ple, Ref. [5] discusses room-temperature magic- angle- 
spinning (MAS) NMR of oxygen, titanium and lead in 
PbZri^^Ti^Oa (PZT, < a; < 1) and the effect of the 
composition (x) on NMR spectra. The present authors 
are unaware of any studies of the longitudinal and trans- 
verse relaxation times of ^°^Pb in PT and PZT below 
room temperature prior to the present work. 



B. Studies of relaxation mechanisms 

Of interest to us here is primarily longitudinal relax- 
ation of nuclear spins in compounds where all electron 
angular momenta are paired (with a possible exception of 
paramagnetic impurities). Several mechanisms can lead 
to such relaxation [101 E] • These processes typically oc- 
cur as a result of the modulation of the crystal field or lig- 
and field due to lattice vibrations (phonons) and can be, 
to some degree, experimentally distinguished by the tem- 
perature and magnetic-field dependence of relaxation. 

There are two types of spin-phonon interaction pro- 
cesses: the "direct process" , in which a spin-flip is 
caused by the interaction with a single lattice phonon, 
and a "Raman process" , involving interaction with two 
phonons. The Raman process has been shown in Ref. 
[12] to be the dominant relaxation mechanism for ^"^^Pb 
nuclei in ionic solids such as lead nitrate above their De- 
bye temperature T^i (typically 100 K to 200 K). The 
nuclear spin-relaxation rate due to this process is pro- 
portional to T^ for T> Td (each phonon interaction con- 
tributes a factor of T in this classical regime) . However 
for T<C Td, this process has a T*" dependence [TT]. (This 
can be derived by integrating the Bose-Einstein energy 
densities of the two phonons over their frequencies up 
to the Debye frequency. Since each energy density con- 
tains a factor of frequency cubed, the resulting depen- 
dence is T^.) The relaxation rate for the Raman process 
is independent of the magnetic field [T?. In Ref. [13] . 
longitudinal spin-lattice relaxation of ^''^Pb in polycrys- 
talline samples of PbMo04 and PbCl2 was investigated 
as a function of temperature down to the liquid-nitrogen 
temperature range. A T^ dependence of 1/Ti was ob- 
served and interpreted as second-order Raman process. 
A rather complete discussion of the theory is found in 



Ref. [12] , where not only the T^ scaling with temperature 
and independence of NMR frequency are reproduced, but 
also the magnitude of the relaxation rates is accounted 
for within an order of magnitude (for lead molybdate, 
lead chloride, lead nitrate, thallium nitrate, thallium ni- 
trite, and thallium perchlorate). The model explains the 
relaxation as caused by a Raman process involving the 
interactions between nuclear spins and lattice vibrations 
via a fluctuating spin-rotation magnetic field. 

Unfortunately, more often than not, impurities are 
responsible for much of the observed nuclear spin- 
relaxation. Back in the early days of NMR, Bloembergen 
fl4J investigated nuclear spin-relaxation of several (non- 
Pb) compounds between 300 and 1 K, and found impuri- 
ties to be the dominant relaxation agent. Spin diffusion 
is essential in affecting relaxation for nuclei that arc far 
away from an impurity. Dependences of the relaxation 
rates ocT and cxT^ were observed at low temperature. 
A further investigation of nuclear spin-relaxation caused 
by paramagnetic impurities was carried out by Blumberg 
|15j . He conducted saturation-recovery experiments and 
found that at short times, there is a \/t dependence of 
the recovery indicative of diffusion. At longer times, the 
recovery is exponential, but the exponent is proportional 
to impurity concentration. 



III. EXPERIMENTAL METHOD 

NMR measurements at the ^'^^Pb frequency were done 
on a 89-mm bore, 11.7 T Oxford magnet using a Var- 
ian NMR Systems (VNMRS) spectrometer console run- 
ning VnmrJ 2. IB software. On this instrument, the pro- 
ton frequency was 500.195 MHz. The RF carrier fre- 
quency was set to 104.481 MHz for PbTiOs (PT) and 
104.495 MHz for PZT. The samples were hghtly packed 
into a 2 cm- long, 0.5 cm-diameter glass tube and excited 
using the solenoidal rf coil of a home-built static probe for 
cryogenic samples. A liquid-helium cryostat surrounded 
the NMR probe and cooled the NMR tube. The tempera- 
ture of the sample was monitored using a resistive sensor 
with a bridge placed in proximity of the glass tube. The 
temporal stability in the temperature control over the 
range 10-200 K was better than 0.1 K; however, tempera- 
ture fluctuations inside the sample were not directly mon- 
itored and could exhibit higher fluctuations and possibly 
some internal temperature gradients. To avoid system- 
atic errors in the relaxation measurements arising from 
possible sample temperature fluctuations, inversion times 
and echo times were ordered randomly for each experi- 
ment and measurements of Ti and T2 were compared to 
those from repeated experiments. We estimate the abso- 
lute error in temperature determination to be less than 
about 1 K. Pulse widths for a 7r/2 rotation ranged be- 
tween 10 and 13 /zs over the temperature range, with the 
dependence mainly due to variations in the coil tuning 
(Q- factor). 

Additional control room-temperature measurements 
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(a) Saturation recovery T1 measurement with spin echo readout 
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FIG. 1: Pulse sequences used for measuring Ti and T2 re- 
laxation times. In both cases, a train of 7r/2-pulses was used 
to saturate the magnetization, followed by a recovery period. 
The Ti measurement used a spin echo readout to avoid coil 
ringing artifacts and several values of Ta were used. T2 mea- 
surements used a Hahn echo for several values of r. The 7r/2 
pulse widths were in the range of 10 — 13 /is. 



with a lead titanate sample were carried out with a 7.1-T 
magnet (proton frequency 300 MHz) with a Varian CMX 
Infinityplus spectrometer, under magic- angle spinning at 
11 kHz. Pulse widths for a 7r/2 rotation on this instru- 
ment was 3 /iS. 

Transverse (T2) relaxation times were measured on the 
11.7-T system using a 7r/2 — t/2 — tt — t/2 — {acquire} 
Hahn echo [Fig. [l] (b)] with echo times (t) at 100 ^is, 
200 /iS, 400 /iS, 800 ^is, 1.6 ms, 4 ms, and 8 ms. The 
preparation consisted of a saturation pulse train followed 
by a 2-min equilibration period for re-polarization of the 
longitudinal magnetization. At a given temperature, this 
provided a consistent initial longitudinal magnetization 
for each point in a T2 measurement. 



IV. SAMPLES 

Experiments were performed on four different sam- 
ples. The PbTiOs powders were obtained from two dif- 
ferent sources: Praxair Specialty Ceramics (chemical pu- 
rity 99.9%, particle size 1 /im), and from the group of 
Prof. M. Jansen at Max Plank Institute for Solid State 
Research in Stuttgart, Germany (chemical purity un- 
known). The powder samples from these two sources 
are denoted as PI and P2, respectively. Measurements 
were also made on a PbTiOs crystal with dimensions 
2.3 mm x 2.8 mm x 1.3 mm. The fourth sample was 
PbZro.52Tio 48O3 powder from Praxair Specialty Ceram- 
ics (chemical purity 99.9%, particle size 1 ^m). 



V. RESULTS AND DISCUSSION 
A. Spin-lattice relcLxation 

Longitudinal (Ti) relaxation times were measured us- 
ing a saturation-recovery sequence shown in Fig. [T| (a). 
The saturation preparation is accomplished by applying 
a train of thirty 7r/2 pulses separated by 30-ms intervals 
used to completely depolarize ^°^Pb spins in the sensitive 
volume despite inhomogeneities in the rf field or inaccu- 
racies in the flip angle. The saturation was checked by 
confirming that no signal was produced in the limit of 
zero recovery time. Signals were acquired with 8 to 15 
recovery times (Tr) logarithmically spaced between ms 
and 5 X Ti. The signals were acquired using a spin echo 
with interpulse delay of 50 ^s to avoid signal contamina- 
tion from probe ringing. 

All the magnetization-recovery data were fit with 
two models: a single exponential M{t) = Mo[l — 
exp(— i/Ti)], and a stretched exponential M{t) = Mq[1 — 
exp(— -y/i/Ti)], where Mq is the equilibrium magnetiza- 
tion, and t is the time elapsed after saturation. The use 
of these models is motivated by the fact that magnetiza- 
tion recovery is well described by a single exponential if 
the dominant spin-lattice relaxation mechanism is a two- 
phonon Raman process [12], or paramagnetic-impurity 
relaxation with fast spin diffusion [16\ If, however, the 
spin-lattice relaxation is due to paramagnetic impurities 
with slow spin-diffusion, then the magnetization recov- 
ers along a stretched exponential curve |16j . It should 
be noted that the chemical shift anisotropy in PbTiOa 
and PbZro.52Tio.48O3 is on the order of 2000 ppm, thus, 
in the applied 11.7-T magnetic field, the local magnetic 
field varies by ~ 200 G. This is much greater than the 
local dipole-dipole interaction field, which is on the order 
of 1 G. This does not mean, however, that spin-diffusion 
is strongly suppressed. The reason for such a broadened 
line is that we are looking at powder. Each powder grain 
has a much narrower line, but, since there is large NMR- 
shift anisotropy due to the chemical shift and diamag- 
netic susceptibility, and the grains are randomly oriented, 
the resultant line is wide. Thus, within each powder crys- 
tallite, spin-diffusion can happen, but the spin-flip cannot 
"hop" between the grains. 

Magnetization recovery data for PbTiOa at T= 10 K 
are shown in Fig. [2] (a). The data points are well fit with 
a single exponential, a stretched exponential gives a poor 
fit. This is the case for all PbTiOa data at all temper- 
atures. The situation is different for PbZro.52Tio.48O3. 
As shown in Fig. |2] (b) , the stretched exponential gives 
a much better fit than a single exponential. This is the 
case for PbZro.52Tio.48O3 at temperatures below 50 K. 
For the data above 50 K, however, single exponential 
provides a better fit. This is an indication that spin- 
lattice relaxation due to paramagnetic impurities, with 
suppressed spin-diffusion, may be the dominant mecha- 
nism for PbZro.52Tio.48O3 below 50 K. 

Longitudinal relaxation rates l/Ti extracted from the 
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FIG. 2: Representative magnetization recovery curves for 
PbTiOa (a), and PbZro.52Tio.4gO3 (b) powders. Two fits are 
shown for each data set: a single exponential and a stretched 
exponential. The FT sample is PI (see text). 



Hahn echo measurement data 
for PbTiOj, P2 sample at 21 K 
- Exponential decay fit witii 
t = (1.62 ± 0.08) ms 
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FIG. 3: The experimentally determined longitudinal relax- 
ation rates vs. temperature. The data were fit to different 
functions above and below 50 K (see text). 



magnetization recovery fits as in Fig. [2] are plotted in 
Fig [31 



B. Spin-spin relaxation 

An example of transverse-relaxation data is shown in 
Fig. |4j Despite of the fact that dipole-dipole interactions 
may lead to non-exponential relaxation [llj . within ex- 
periential uncertainties, all our data were consistent with 
exponential relaxation, so exponential fits were used to 
extract the relaxation times 

Analysis of the transverse-relaxation data at different 
temperatures results in the value T2 = (1.5 ±0.2) ms, for 
all the PbTiOa and the PbZro.52Tio.48O3 powder sam- 
ples. No temperature variation of T2 is observed within 
the experimental uncertainty. (No T2 data were taken 
for the crystal.) 

The measured T2 value is within a factor of two of the 
R:!0.7-ms decay time estimated for PT in Ref. [I] assuming 



FIG. 4: An example of data and exponential-decay fit used 
to extract the T2 values. 



relaxation due to dipole-dipole interaction. The fact that 
the relaxation rate is independent of temperature is also 
consistent with relaxation via this mechanism. 



VI. INTERPRETATION OF Ti RESULTS 
A. High-temperature range 

In the temperature range of 50-293 K, the observed 
spin-lattice relaxation rates for all our samples follow 
the dependence. This behavior is characteristic of a 
two-phonon Raman process [TTl [12l [13l [17] , a relaxation 
mechanism investigated and explained in terms of modu- 
lations in the spin-rotation interaction by Dybowski and 
co-workers [T2J [T^ . Our results, which are the first mea- 
surements below room temperature on PT and PZT per- 
ovskite ferroelectrics, confirm that this relaxation mech- 
anism is, indeed, the dominant process over the 50-293 K 
temperature range. We note that Debye temperatures of 
ionic crystals are in the range 150-350 K, whereas our 
results exhibit dependence down to 50 K. 



B. Low-temperature range 

At low temperatures (below 50 K) the dependence 
is no longer a good fit to the data shown in Fig. |3] An- 
other mechanism dominates spin-lattice relaxation in this 
regime - this mechanism is most likely the relaxation due 
to paramagnetic impurities in the samples. The ■^'^^Pb 
nuclear magnetic sublevels are split by an energy corre- 
sponding to about 5 mK in the applied field of 11.7 T. 
Phonons which correspond to that same frequency would 
normally cause ^''^Pb nuclear spin-flips via the direct pro- 
cess. However, since 5 mK is much smaller than the 
sample temperature, the energy density of such phonons 
is greatly suppressed. An impurity with uncompensated 
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electron spin, however, has a much larger magnetic mo- 
ment. In fact, its magnetic sublevels are split by about 
8 K in the applied magnetic field, which is on the order 
of our samples temperature (at the lowest temperatures 
where the data were taken). Therefore such an impu- 
rity can interact with the lattice via the direct phonon 
process, i.e. absorption or emission of a single phonon 
at the Larmor frequency of the electronic spin ws- For 
a Kramers doublet, the relaxation rate of the electronic 
spin is given by the Kramers relation [17]: 
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where v is the velocity of sound in the crystal, p is the 
density and ^ is a dimensionless coefiicient that depends 
on the structure of the impurity. For non-Kramers sites, 
there is a different expression for the relaxation rate, how- 
ever, the temperature dependence of T^^ is the same as 
in Eq. Q 110^. 

Spin-lattice relaxation of nuclei (I) in the vicinity of an 
impurity can be described by a dipolar coupling S • I? • I, 
where T) is, & tensor describing the magnetic dipolar cou- 
pling between the electron and the nucleus (separated by 
a distance r), while S - I? can be viewed as a random fluc- 
tuating local field produced by the electron and seen by 
the nucleus. At low temperatures where motion is frozen, 
the fluctuations are mainly due to the spin flipping of the 
electron, the timescale of which is the relaxation time Ti^. 

When the power spectrum of these fluctuations con- 
tains non-vanishing values at the Larmor frequency ui 
of the nucleus, nuclear relaxation will take place. An ar- 
gument based on the freezing of relative motions at low 
temperatures [T7] leads to the following estimate for the 
nuclear spin-lattice relaxation rate of I due to fluctua- 
tions in S: 
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where Pg — ta,'Dh{hwg /2kBT) is the electronic thermal 
equilibrium polarization and Be = fif,Sr~^ gives the or- 
der of magnitude of the local magnetic field seen by the 
nucleus due to electron magnetic moment /ig. A further 
level of refinement involves integration of the spatially- 
dependent rates Tin~^{r) in close proximity of the im- 
purities |17j . Finally, mutual spin flip-fiops of the nu- 
clei cause this relaxation mechanism to spread its effects 
throughout the lattice. The (temperature-independent) 
rate of spin diffusion and density of impurities determine 
whether relaxation is best described by monoexponential 
or stretched-exponential behavior 16]. 

The low-temperature (< 50 K) spin-lattice relaxation 
data were flt with the function 
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which is just Eqs.Qand Q combined. The results are 
shown with curved lines in Fig. [3] The fit contained one 



free parameter, A. The value hujs/2kB = 8 K was calcu- 
lated from the applied magnetic field of 11.7 T (assum- 
ing the impurity g-factor of 2), and fixed for the fits to 
PbTiOs PI powder data and the PbZro.52Tio.48O3 data. 
It was found, however, that the PbTiOs P2 powder data 
were fit much better with the value hujs/2kB — 16 K, 
i.e. the paramagnetic impurities in this sample have the 
magnetic moment of 4 Bohr magnetons. We can con- 
clude that the model described above gives a satisfac- 
tory description of the data below 50 K. The PI PbTiOa 
powder seems to have the smallest concentration of para- 
magnetic impurities, while the PbZro.52Tio.48O3 and the 
P2 PbTiOa powders have larger impurity concentrations, 
and the dominant impurity in PbTiOa P2 powder has the 
g-factor of 4. 



C. Non-exponential relaxation 

We note that although wc fit the relaxation data for 
PZT to a stretched exponential for the reasons discussed 
above, equally good fits can be obtained by using a biex- 
ponential, where the rapid component is roughly half the 
amplitude of the slower component. Possible additional 
mechanisms which could lead to the non-exponential 
behavior include the following. 1) The PZT subunits 
PbZrOa and PbTiOa may exhibit different spin-rotation 
interactions that modulate the relaxation rate Tj~^ of 
the Raman process Pil]. 2) The stretched exponential 
may arise from a distribution of relaxation times, due to 
some degree of anisotropy in the relaxation rate, how- 
ever, note that such relaxation behavior is not seen in 
PT. Unfortunately, the relaxation data alone arc insuf- 
ficient to determine the source of the deviation from a 
single exponential. 



VII. IMPLICATIONS FOR THE 
SCHIFF-MOMENT EXPERIMENT 

The results shown in Fig. |3]are encouraging for the on- 
going development of the experiment to search for the P- 
and T-violating Schiff moment of ^°^Pb in ferroelectrics. 
As expected, the spin-lattice relaxation time increases 
dramatically as temperature is lowered. The presence 
of paramagnetic impurities leads to approximately lin- 
ear temperature dependence of the relaxation rate below 
50 K: 1/Ti (xT, which is much slower than the T^ depen- 
dence, expected if the Raman phonon relaxation mech- 
anism were to dominate. This may prove advantageous 
for the design of the Schiff-moment experimental search, 
which relies on achieving a high degree of nuclear spin po- 
larization in high magnetic field and at low temperature, 
prior to nuclear demagnetization to reach ultra-low spin 
temperature [2j. Indeed, with Ti on the order of an hour 
at lattice temperature of 10 K, such polarization can be 
established on this time scale, which is not prohibitively 
long. After demagnetization, the spin-temperature stays 
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low for time also on the order of Ti ^ hour, which 
is when, for example, the SchifF moment-induced mag- 
netization can be measured by precision magnetometry 
methods. 

Two main questions remain to be answered in future 
investigations: what is the spin-lattice relaxation time 
in low magnetic fields and low spin temperatures, and 
whether or not the presence of paramagnetic impurities 
introduces serious systcmatics into the SchifF moment- 
induced magnetization search. We are currently devel- 
oping experiments to address both of these questions. 

VIII. CONCLUSION 

In conclusion, we have presented the first experimental 
study of relax;ation properties of ^'^^Pb in PT and PZT 
below room temperature. We find that above T» 50 K, 
longitudinal relaxation rate follows the dependence 
characteristic of the two-phonon Raman process. On 
the other hand, as the temperature is decreased below 
Tw 50 K, the longitudinal relaxation rates drop slower 
than ocT^ (as opposed to ocT^ expected for the Raman 
process), and the relaxation is probably due to a di- 
rect process associated with paramagnetic impurities and 
nuclear-spin diffusion. While the longitudinal relaxation 
times Ti vary between several seconds and over an hour in 
the temperature range between 290 and 10 K, the trans- 



verse relaxation time T2 is found to be « 1.5 ms for all 
temperatures and all powder samples studied. The ob- 
tained results provide an important input in the design 
of the experiments to search for P,T-violating effects in 
solid ferroelectrics. 
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